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PREFACE 

The  Acoustics  and  Seismics  Laboratory  of  The  University  of  Michigan's  Insti¬ 
tute  of  Science  and  Technology  has  conducted  research  in  seismology  for  several 
years  under  the  sponsorship  of  various  agencies  such  as  the  U.  S.  Air  Force  Office 
of  Scientific  Research,  Air  Force  Technical  Applications  Center,  Air  Force  Cam¬ 
bridge  Research  Laboratories,  and  National  Science  Foundation.  Ehiring  the  course 
of  this  extended  program  the  laboratory  has  obtained  a  library  of  seismic  data, 
which  is  available  for  study. 

This  report  investigates  an  auditory  analysis  method  applied  to  some  of  the 
data,  and  summarizes  2  1/2  years  of  study,  ending  15  February  1964,  sponsored  by 
the  Air  Force  Office  of  Scientific  Research  under  Contract  AF  49(638)- 1079  as  purt 
of  the  Advanced  Research  Projects  Agency’s  VELA  UNIFORM  program. 

Field  measurements  were  obtained  primarily  under  Air  Force  Contracts 
AF  49(638)- 1 1 70,  administered  through  the  Air  Force  Office  of  Scientific  Research; 
AF  19(6401-8809  and  AF  19(6041-6642,  administered  through  the  Air  Force  Cam¬ 
bridge  Research  Laboratories;  and  AF  49(6381-911,  administered  through  the  Air 
Force  Technical  Applications  Center. 

This  report  was  prepared  by  G.  E.  Frantti  as  the  contractually  required  finai 
report  on  Contract  AF  49(6381-1079  and  as  a  technical  report  on  Contract  AF  49(6381- 
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INVESTIGATION  OF  AUDITORY  DISCRIMINATION  OF 
SEISMIC  SIGNALS  FROM  EARTHQUAKES  AND  EXPLOSIONS 


ABSTRACT 

Magnetic  tape  recordings  of  short-period  seismic  signals  from  approximately 
200  earthquakes  and  explosions  were  time-compressed  by  a  factor  of  up  to  512  to 
shift  seismic  frequencies  to  the  audible  range.  These  seismic  data  Include  the 
Inhomogeneities  Introduced  by  substantial  variations  In  the  locations  of  sources 
and  receivers  (world-wide),  propagation  path  length  (32  to  7000  km),  and  source 
magnitude  (M  '  0.5  to  M  =  6.5).  Subjects  were  trained  with  a  representative  set 
of  the  "seismic  sounds."  Auditory  experiments  were  conducted  to  determine  the 
ability  of  the  human  auditory  system  to  distinguish  between  seismic  signals  from 
earthquakes  and  explosions.  The  results  of  the  experiments  suggest  that  a  trained 
listener  can  Identify  approximately  two  thirds  of  the  seismic  sounds  presented. 


1 

INTRODUCTION 

Soismologicai  research  has  been  conducted  for  several  years  by  the  Selsmlcs  Group  at  the 
Acoustics  and  Selsmlcs  Laboratory  of  The  University  of  Michigan’s  Institute  of  Science  and 
Technology  During  this  extended  program  the  laboratory  hits  accumulated  a  library  of  mag¬ 
netic  tape  recordings  of  short-period  seismic  signals  (0.5  to  500  cps)  from  such  sources  as 
earthquakes,  nuclear  detonations,  high-explosive  tests,  and  quarry  blasts. 

One  objective  of  the  research  is  to  Investigate  techniques  which  might  aid  In  distinguishing 
between  earthquake  seismograms  and  explosion  seismograms.  This  problem  Is  of  primary 
concern  in  the  VELA  UNIFORM  program,  and  has  widespread  Interest  within  the  soismologicai 
discipline.  Important  practical  progress  has  been  made  by  the  application  of  many  analytical 
techniques  to  seismograms,  but  It  Is  apparent  that  new  approaches  are  still  needed. 

The  experiments  described  in  this  report  were  conducted  to  explore  the  ability  of  the  human 
ear  lo  discriminate  between  earthquakes  and  explosions  by  their  audio  signature.  Time  com¬ 
pression  of  the  tape-recorded  signals  is  used  lo  shift  seismic  frequencies  Into  the  audible  range. 
Should  such  a  subjective  technique  be  found  useful.  It  would  represent  a  rapid  means  of  mon¬ 
itoring  seismometer  outputs. 

After  conducting  auditory  experiments  at  the  Bell  Telephone  Laboratories,  Speeth  con¬ 
cluded  that  listeners  identified  two  classes  of  events  explosions  and  earthquakes,  In  over  90% 
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of  the  cases  presented  them  (1).  However,  experimental  data  available  for  that  study  were 
severely  limited.  In  later  experiments  [2)  the  level  of  identification  was  less  than  90  f,,  but  the 
results  were  still  not  definitive.  Thus,  though  the  application  of  the  auditory  technique  to  seis¬ 
mic  data  is  not  new,  it  has  not  been  adequately  explored. 

The  auditory  research  program  reported  here  proceeded  in  three  phases: 

Phase  I  -  Equipment  was  purchased,  built,  and  integrated  to  develop  an  instrumentation 
system  for  time-compressing  the  tape-recorded  seismic  signals. 

Phase  II  -  Seismic  recordings  were  processed  through  the  tape  speed-up  system,  and  an 
auditory  data  library  was  built  up. 

Phase  III  -  The  auditory  library  was  used  to  train  subjects  and  carry  out  a  set  of  auditory 
experiments. 


2 

DISCUSSION  AND  RESULTS 

2.1.  SEISMIC  DATA 

About  200  seismograms  of  known  seismic  events  were  used  In  this  study:  examples  are 
shown  in  Figure  1.  They  were  equally  distributed  iietwcen  earthquakes  and  explosions.  Re¬ 
cordings  were  made  in  the  ficid  with  Ampex  FM  magnetic  tape  recorders,  operated  at  3  3/4  Ips 
and  using  seismic  amplifiers  built  at  The  University  of  Michigan  (3-db-down  points  0.5  and 
800  cps).  At  most  stations  the  ficid  detector  was  a  Hall-Sears  (US- 10)  2-cycle  geophone,  but 
some  of  the  data  were  obtained  with  I -second  Hen l off  and  l-sccond  Wlllmore  seismometers. 
Typical  ficid  equipment  are  shown  in  Figure  2.  In  a  relatively  recent  Innovation,  a  capacitor 
across  tne  IIS-10  is  used  to  decrease  the  sensitivity  to  high-frequency  noise.  Three  earthquake 
recordings  made  with  short-period  instruments  were  obtained  from  the  Geotechnical  Corpor¬ 
ation. 

The  set  of  fieid  data,  then,  was  obtained  with  an  approximately  consistent  recording  system. 
Seismometers  were  generally  set  up  on  licdrock  and  shallow  burled  to  reduce  local  effects. 

Local  and  near-regional  earthquakes  were  recorded  in  a  number  of  locations  in  North 
America,  as  well  as  in  Chile,  Hawaii,  Crete,  American  Samoa,  and  Puerto  Hlco.  Also,  earth¬ 
quakes  with  epicenters  in  Mexico,  the  Kurlic  Islands,  Peru,  the  Azores,  and  the  Aleutians  were 
recorded  at  teleseismic  distances.  The  propagation  distance  ranged  from  32  to  7000  km.  Epi¬ 
centers  were  in  both  water-covered  and  continental  areas.  The  data  cover  a  range  in  source 
magnitude  of  M  =  0.5  to  M  =  6.5. 

The  explosions,  ail  recorded  in  North  America,  include  nuclear  detonations,  high-explosive 
tests,  underwater  shots,  and  quarry  blasts.  Source  size  varied  from  one-half  ton  to  about  100 
kliotons  (equivalent  in  TNT),  and  the  propagation  distance  varied  between  50  and  3100  km. 
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FIGURE  1.  EXAMPLES  OK  SKISMOC.RAMS  USED  IN  THIS  STUDY  Wideband  awl  passband  traces  arc  shown  Playback  gain  is 

Indicated  In  db  The  seismometer  was  a  two-cycle  Hall-Scars 
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This  data  population  Is  sucii  tKat  II  Is  possible  to  assemble  subsets  of  events  with  (I)  dis¬ 
tance  held  constant  and  event  magnitude  varied,  (2)  distance  the  variable  and  magnitude  held 
approximately  constant,  or  (3)  complete  data  Inhomogcncity  with  all  pirameters  varying  ran¬ 
domly. 

2  2.  DATA  PROCESSING 

The  seismic  recordings  used,  which  are  described  in  detail  In  the  appendix,  arc  on  mag¬ 
netic  tape,  a  form  convenient  for  time-compressing  the  data.  Figure  3  is  a  diagram  of  this 
speed-up  process,  and  Figure  4  shows  the  main  components  of  the  system.  Original  field  tapes 
are  initially  played  through  a  standard  data  top  plate  to  double  the  tape  speed.  The  tape  is  then 
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demodulated  and  passband  filtered  at  0.2  to  40  cps  (0.1  to  20  cps  in  real  time)  to  eliminate  some 
of  the  noise.  The  next  stage  consists  of  a  d-c  line  amplifier  which  normalizes  the  data;  that  is, 
the  maximum  amplitude  on  a  seismogram  is  adjusted  to  one  volt  rms.  The  normalized  signals 
are  then  recorded  at  3  3/4  ips  on  track  one  of  a  2-track  1/4-inch  loop  machine  (loop  length  75 
feet)  built  at  The  University  of  Michigan.  The  second  track  of  the  1/ 4-inch  ioop  contains  a  con¬ 
trol  signal  which,  at  the  time  the  tape  splice  approaches  the  reproduce  head,  switches  to  a  sec¬ 
ond  reproduce  head  advanced  on  the  ioop,  thus  eliminating  the  transient  otherwise  introduced 
by  the  spiicc.  Output  from  the  1/4-tnch  machine  is  played  at  15  ips,  so  that  at  this  point  in  the 
processing  an  8x  speedup  is  achieved.  The  tape  Is  again  filtered,  at  0.2  to  160  cps,  and  then 
recorded  at  1  7/8  ips  on  one  track  of  a  one-inch,  14-track  Ampex  FL-200  B  Loop  Machine 
(ioop  tength  150  feet).  The  output  of  this  machine  is  operated  at  15,  30,  or  60  ips  to  give  over¬ 
all  time  compressions  of  64,  128,  or  256,  respectively.  These  time  compressions  can  be  in¬ 
creased  by  generating  higher  speedups  in  eariler  stages  of  the  process  for  example,  in  the 
1/4-lnch  loop  machine.  The  output  of  the  one-inch  ioop  is  preampitfied,  passed  through  a  pair 
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FIGl'RE  I.  PRIMARY  COMPONENTS  OE  THE  TIME  COMPRESSION  SYSTEM  f  mm  the  left. 
I  '4 -Inch  loop  machine,  I -inch  loop  machine  (door  closed),  standard  data  top  plate,  amplifiers 
electronics,  and  power  supplies. 
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of  Mar-anil  70-walt  power  amplifiers,  and  finally  displayed  through  a  pair  of  AR-3  speakers  In 
an  acoustically  favorable  room.  A  second  output  on  the  preamplifier  stage  is  recorded  on  f/4- 
inch  reels  for  storage  and  preparation  of  training  tapes. 

The  final  product  of  the  time-compression  process,  then,  is  a  set  of  one-inch,  M-lrack, 
f  50-foot  tape  loops.  Each  loop  contains  M  seismic  events  (one  per  channel),  and  each  event  is 
repeated  f  6  times  on  Its  channel.  A  head  switching  panel  associated  with  the  one-inch  loop 
machine  enables  rapid  selection  of  the  desired  channel. 

2.3.  BACKGROUND 

The  signals  recorded  on  one-inch  loops  convey  complex  audio  Information  to  the  listener. 
This  Information  consists  of  multiple  stimuli  which  arise  from  interrelationships  among  vari¬ 
ous  phases  of  seismic  waves.  Each  stimulus  must  be  statistically  related  to  a  particular  subset 
of  events  (earthquakes,  EQ,  or  explosions,  EX).  Signal  identification  is  then  a  case  of  testing 
statistical  hypotheses,  where  the  listener  adopts  some  system  of  optimizing  to  select  the  par¬ 
ticular  subset  a  sound  belongs  to. 

For  problems  In  signal  detection  it  has  been  shown  (see.  for  example,  Reference  3)  that 
the  best  receiver  (listener)  Is  the  one  which  calculates  a  likelihood  ratio  for  each  Input,  fn  this 
way  the  events  are  mapped  onto  a  single  axis.  Then  the  listener  can  decide  which  subset  an 
input  was  taken  from  by  applying  some  rule  pertaining  to  likelihood,  or  to  some  point  on  the 
axis,  ff  there  arc  many  stimuli  associated  with  a  given  Input,  as  In  the  seismic  case,  It  Is  nec¬ 
essary  to  calculate  a  separate  likelihood  ratio  for  each  criterion.  The  possibility  of  there  being 
as  many  criteria  as  there  are  listeners  is  perhaps  not  remote. 

Since,  In  fact,  It  Is  not  possible  to  tdonllfy  all  the  stimuli,  it  is  difficult  to  conduct  experi¬ 
ments  based  on  an  optimizing  procedure  which  comprehensively  investigate  the  auditory  dis¬ 
crimination  capability  in  the  seismic  case.  However,  we  anticipate  that  the  combined  effect  of 
multiple  criteria  can  be  readily  explored  by  evaluating  the  l>ehavtor  of  listeners  in  a  'sound 
identification"  program  according  to  the  percentage  of  correct  decisions. 

Tanner  [<1  expands  the  theory  of  auditory  detection  of  signals  to  treat  a  simple  signal  rec¬ 
ognition  problem,  with  the  support  of  auditory  experiments.  His  argument  will  not  be  repro¬ 
duced  here,  but  it  does  demonstrate  that  the  theory  of  statistical  decision  is  applicable  to  the 
problem  of  recognition. 

The  available  data  In  this  study,  then,  consist  of  two  subsets:  population  EQ  (earthquakes) 
and  population  EX  (explosions).  A  sample  Is  drawn  from  one  of  these  two  distributions,  and  the 
observer  must  Identify  the  distribution  from  which  it  was  taken.  We  define  a  decision  region, 
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A,  such  that  II  an  observation  tails  In  this  region  the  listener  accepts  the  hypothesis  that  the 
event  Is  In  the  set  EQ.  It  the  observation  lalls  In  the  region  CA,  the  listener  accepts  the  hy¬ 
pothesis  that  the  event  Is  In  the  set  EX.  Thus  we  have  tor  the  possible  ou.comcs  ot  the  experi¬ 
ment 

P(FQ  •  A)  •  P(EQ  •  CA)  •  P(EX  •  A)  ♦  P(EX  •  CA)  »  1  (1) 

where  the  terms  on  the  left,  In  order,  are  the  probabilities  that  (a)  given  EQ  the  listener  will 
accept  It  as  EQ.  (b)  given  EQ  he  will  accept  it  as  EX,  (c)  given  EX  he  will  accept  It  as  EQ,  and 
(d)  given  EX  lie  will  accept  it  as  EX.  Now  the  conditional  probabilities  associated  with  the  ex¬ 
periment  are  such  that  the  following  statements  may  also  be  made: 


P(EQ)  •  P(EX)  *  1 

(2) 

P(EQ)  *  P(EX)  -  0.5 

(3) 

peq(A>  *  peq(ca)  - 1 

W) 

pex(A)  •  pex(ca)  * 1 

(5) 

Since  Equations  4  and  5  equate  to  one,  we  may  describe  the  perlormance  of  listeners  with  Just 
two  probabilities,  for  example  PEQ(A)  and  PEX(A).  The  lormcr  Is  the  probability  that  a  lis¬ 
tener  will  Identity  an  EQ  correctly.  The  latter  Is  the  probability  that  he  will  Identity  an  EX  as 
an  EQ.  In  the  literature  these  are  referred  to  as  hit  rate  and  false  alarm  rate,  respectively, 
and  by  standard  dcllnation  In  lormal  probability  theory  they  are 


WAI  ■  J 

>(x,dx 

■  J 

['ex00*1 

(6) 


where  tho  Integration  is  over  all  points  In  region  A,  and  the  integrands  arc  the  probability  den¬ 
sity  lunctlons  of  sample  variable  x  for  the  cases  when  x  Is  drawn  lrom  populations  EQ  and  LX, 
respectively. 

In  terms  ol  the  density  functions  above,  likelihood  ratio  can  be  expressed  as 


f(x) 


'eq(x) 

'ex00 


(7) 


Suppose  we  ask  the  listener  to  maximize  his  hit  rate  relative  to  his  lalse  alarm  rate,  that  is, 


PEq(A)  -  w  Pex(A)  =  max 


(8) 
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which  Is  the  same  as  maximizing  the  Integral 

/  lfEQ(x)  '  ’**’  'eX(x)'  ^  ‘  maX  (9) 

Then  we  see  that  the  decision  rule  Involves  a  ratio  criterion.  Specifically  this  says  choose 
region  A  such  that 

f(x)>*-  no) 

being  the  weighting  factor  introduced  In  Equation  8. 

Consider  the  probability  density  functions  in  the  seismic  recognition  problem  to  be  as 
shown  In  Figure  5.  The  extent  to  which  seismic  signals  can  be  discriminated  Is  depleted  by 
the  separation  of  the  means  of  distributions  EX  and  EQ.  If  the  listener  chooses  a  criterion  cut 
at  some  point  C  (related  to  a  specific  value  of  likelihood  ratio),  then  the  cross-hatched  area  Is 
a  measure  of  hit  rate  and  the  stippled  area  Is  a  measure  of  false  alarm  rate.  As  the  cut  Is 
moved  to  different  positions  on  the  decision  axis,  the  conditional  probabilities  change  and  a  set 
of  points  relating  hit  rate  and  false  alarm  rate  Is  generated.  The  relationship  between  such 
points  Is  known  as  an  ROC  (receiver  operating  characteristic)  curve.  If  these  points  plot  In  a 
straight  line  with  slope  one  when  presented  on  double- probability  paper,  the  probability  den¬ 
sities  are  normally  distributed  and  have  equal  variance. 

The  value  of  an  ROC  presentation  Is  that  It  completely  describes  an  observer’s  performance 
In  "sensory**  experiments  where  decisions  are  based  on  ratio  criteria,  ftecause  of  the  Inherent 


SUBJECT'S  CRITERION 


FIGURE  5.  HYPOTHETICAL  PROBABILITY  DENSITY  FUNCTIONS  OF  EQ  AND  FX  AS  A 
FUNCTION  OF  LISTENER  CRITERION.  An  arbitrary  cut  C  related  to  likelihood  ratio  Is 
shown.  The  cross-hatched  area  represents  P|tq(A).  and  the  stippled  area  represents  Pyx*''* 
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complexity  o(  seismic  signals,  we  choose  to  measure  the  auditory  performance  on  the  basis  of 
percentage  of  correct  decisions,  always  with  the  underlying  objective  of  obtaining,  if  possible, 
likelihood  ratio  information  of  the  type  described.  The  experiments  are  described  In  the  fol¬ 
lowing  section. 

2.4.  EXPERIMENTAL  PROCEDURE 

The  auditory  program  consisted  of  two  experiments,  Series  A  and  Series  B.  Series  A  w-as 
carried  out  at  the  time  when  approximately  125  seismic  events  were  time-compressed.  This 
series,  which  lasted  two  months  included  training  sessions  and  eigiil  test  sessions.  It  was  con¬ 
ducted  (1)  to  obtain  an  early  indication  of  results  (after  preliminary  training);  (2)  to  investigate 
such  factors  as  auditory  stimuli,  speedup,  filtering,  and  listening  to  tapes  played  backwards; 
and  (3)  to  increase  the  listener's  skill.  Series  B  was  carried  out  after  the  entire  sample  (about 
200  events)  was  time-compressed.  This  series,  which  also  lasted  two  months,  included  training 
sessions  and  16  test  sessions.  It  was  conducted  (1)  to  obtain  definitive  results  of  the  percentage 
ol  correct  decisions  that  a  trained  listener  couid  achieve  in  an  identification  program  of  random, 
inhomogeneous  seismic  sounds;  (2)  to  generate  data.  If  possible,  to  describe  the  observer’s 
performance  by  an  ROC  curve;  and  (3)  to  compare  auditory  results  with  analytical  results  ob¬ 
tained  with  the  same  raw  data.  Listeners  were  trained  both  before  and  during  both  series  of 
experiments.  For  tills  purpose  data  other  than  the  lest  data  were  used. 

In  all  experiments  the  subjects  (from  9  to  13  per  session)  were  sealed  In  an  acoustically 
lavorahic  room  In  which  two  Alt- 3  speakers  w-ere  located.  All  other  related  equipment  re¬ 
mained  in  a  separate  room.  Eacli  event  played  through  the  speakers  was  preceded  and  followed 
by  a  short  segment  ol  tiackground  noise  w  hich  conditioned  the  ear  to  the  Impending  signals. 

Each  test  sound  was  presented  lour  times,  the  listeners  having  decided  that  four  was  optimum. 
Although  the  a  priori  protubliity  of  a  given  event’s  being  an  EQ  or  EX  was  0.5,  the  order  of 
presentation  of  unknown  sounds  was  random  (based  on  random  numlicr  tables).  In  the  various 
test  sessions,  data  were  grouped  to  give  the  following  combinations  of  parameters:  (1)  distance 
the  constant  and  event  magnitude  a  variable,  (2)  distance  the  variable  and  event  magnitude  ap¬ 
proximately  constant,  and  (3)  complete  data  Inhomogonelty  with  alt  parameters  varying. 

2.4.1.  SERIES  A.  Series  A  was  a  simple  two-choice  experiment  In  which  listeners  re¬ 
sponded  either  EQ  or  EX.  A  total  ol  21  observers  took  part  In  this  experiment,  although  no  more 
than  13  were  present  in  any  given  test  session.  Subjects  for  both  series  were  staff  members  of 
the  Acoustics  and  Seismics  Laboratory,  including  engineers,  technicians,  secretaries,  and  stu¬ 
dent  assistants. 
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Early  phases  of  this  experiment  Investigated  the  effect  of  filtering.  Seismograms  were 
exposed  to  different  passband  filters,  and  results  were  assessed  in  the  listening  program.  For 
most  of  the  recordings,  a  decrease  in  the  high-frequency  cutoff  to  10  cps  did  not  result  in  any 
perceptual  change  in  the  audio  output.  High-frequency  seismograms  from  local  earthquakes 
recorded  in  Chile,  and  certain  other  recordings  of  near  events  (<  150  km),  were  the  main  ex¬ 
ceptions.  In  fact,  a  high-frequency  limit  of  5  cps  appeared  to  be  satisfactory  for  most  sounds 
obtained  at  distances  exceeding  about  300  km.  The  low-frequency  cutoff  could  not  be  raised 
above  0.5  cps  without  introducing  noticeable  loss  of  information  in  the  output.  The  system  re¬ 
sponse  for  data  gathered  for  this  study  decreases  rapidly  below  0.5  cps.  We  infer  from  these 
observations  that  the  critical  frequency  window  for  a  majority  of  real-time  seismic  data  is  ap¬ 
proximately  0.5  to  5  cps.  To  allow  a  margin  of  safety  in  this  study  we  chose  to  filter  the  orig¬ 
inal  seismograms  at  0.1  to  20  cps,  which  is  limited  enough  to  exclude  much  of  the  spurious 
noise.  In  exceptional  cases  the  high-frequency  cutoff  was  appropriately  increased  above  20  cps 
to  include  aii  short-period  seismic  energy. 

A  trial  and  test  procedure  similar  to  that  described  above  was  conducted  to  explore  the 
effect  of  different  time-compression  factors,  which  ranged  from  one  (real  time)  to  512.  Chi  the 
basis  of  listeners'  performance,  it  was  determined  that  a  time-compression  factor  of  128  is 
optimum  for  events  recorded  in  the  approximate  range  from  200  to  1000  km.  This  result  sug¬ 
gests  that  in  the  speed-up  process  a  minimum  of  about  1/4  second  should  be  maintained  iietween 
the  predominant  P  and  predominant  shear-surface  arrivals,  and  that  the  overall  signal  duration 
should  lie  less  than  about  2  seconds. 

A  resume  of  results  for  eight  test  sessions  in  Series  A  is  shown  in  Table  I.  P(c)  is  the 
mean  percentage  of  correct  identifications  for  the  group  of  listeners  present  at  each  test  session, 
and  o  is  the  standard  deviation  of  the  mean.  The  values  of  P(c)  for  tape  loops  number  4  and  11 
in  this  table  are  examples  of  the  type  of  evidence  which  supports  the  time-compression  vs. 
distance  relationship.  Seismic  events  on  those  loops  are  in  the  ranges  214  to  284  km  and  270 
to  780  km,  respectively,  for  which  the  factor  128  is  optimum. 

Three  of  the  data  loops  were  run  backward  in  some  tests,  with  the  idea  that  the  late-arriving 
coda  of  waves  on  a  seismogram  might  be  a  preferred  "first  arrival  to  the  ear  over  the  much 
more  abrupt  first  P  arrival.  The  results  for  backward  listening,  tabulated  with  a  symbol  (13)  in 
the  P(c)  coiumn  of  Table  I,  indicate  that  backward  listening  does  not  Improve  discrimination. 
However,  it  turns  out  that  subjects  become  confused  when  forced  to  listen  to  seismograms  in 
both  directions.  For  example,  one  of  the  decision  criteria  is  "attack  time,  which  is  longer  for 
earthquakes  intheforwarddirectlon  but  longerfor  explosions  in  the  opposite  direction.  Thus,  to 
properly  exploit  reversed  piaylng  of  seismic  sounds,  one  should  train  a  separate  group  of  subjects 
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TABLE  1.  RESULTS  FOR  SERIES  A 
No.  ol 


Session 

Loop 

Observers 

Responses 

Speedup 

P(c) 

a 

, 

5 

9 

504 

128X 

0.57 

0.08 

2 

4 

12 

864 

64X 

0.60 

0.08 

128X 

0.68 

0.07 

256X 

0.59 

0.08 

3 

8 

13 

910 

128X 

0.58 

0.08 

4 

10 

13 

910 

128X 

0.65 

0.14 

5 

10 

11 

770 

128X 

0.62(B) 

0.08 

6 

6 

10 

400 

64X 

0.58 

0.10 

64X 

0.55(B) 

0.10 

7 

11 

400 

64X 

0.63 

0.07 

64X 

0.61(B) 

0.08 

7 

9 

11 

528 

64X 

0.82 

0.09 

11 

11 

528 

128X 

0.60 

0.12 

8 

11 

12 

432 

256X 

0.51 

0.06 

From  the  data  In  Table  I,  the  mean  P(c)  lor  Series  A  Is 

P(c)  •  0.64  *0.07  (II) 

This  average  Is  based  on  all  data  exclusive  ol  (a)  results  lor  non-optimum  speedups  In  the  case 
ol  loops  4  and  II  and  (b)  scores  labeled  (B). 


2.4.2.  SERIES  B.  Scries  11  was  a  torccd-cholcc  experiment  In  which  listeners  responded 
In  one  ol  lour  categories  and  attempted  to  maximize  the  value  ol  their  response.  The  response 
categories  and  their  ratings  are  given  In  Table  H.  It  Is  considered  that  II  an  observer  Is  re¬ 
quired  to  make  a  discrete  response,  as  In  Scries  A,  then  some  Inlormatlon  available  to  him  Is 
discarded  In  the  decision  process.  In  using  a  rating  system,  the  listener  Is,  In  a  sense,  making 
a  statement  ol  likelihood,  and  thus  more  of  the  auditory  Information  Is  retained  In  his  response. 
This  rating  system  was  chosen  lor  Its  simplicity. 


TABLE  11.  VALUES  AND  COSTS  ASSOCIATED  WITH 
FOUR  RESPONSE  CATEGORIES 


Rating 


Response 

Correct  Decision 

Incorrect  Decision 

I'm  sure  It  is  EQ 

+2 

-2 

1  think  it  is  EQ 

+  1 

-1 

1  think  it  is  EX 

+1 

-1 

I'm  sure  It  is  EX 

+2 

-2 
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A  total  of  19  observers  took  part  In  the  Series  B  experiment.  No  more  than  12  were  present 
in  any  given  session,  and  12  of  the  19  also  took  part  in  Series  A. 

Table  1(1  outlines  11  sessions,  in  which  16  data  loops  were  tested.  In  this  table  P’(c)  refers 
to  the  weighted  average  of  the  group  percentage  of  correct  Identifications  for  each  data  loop. 

In  Table  111  the  results  of  all  19  subjects  were  used.  Decisions  In  the  rating  category  *2  were 
given  a  weighting  of  1,  and  decisions  In  the  rating  category  *1  were  given  a  weighting  of  0.75 
The  mean  of  these  group  averages  Is 

P‘(c)  =  0.618  tO.054  (12) 

Comparison  of  equations  11  and  12  indicates  that  indeed  there  Is  no  significant  difference  be¬ 
tween  the  results  for  Series  A  and  Series  B. 


TABLE  Ill.  RESULTS  FOR  SERIES  B 
Number  of 


Session 

Loop 

Observers 

Responses 

Speedup 

P’(c) 

1 

9 

8 

336 

64* 

0.651 

2 

7 

12 

504 

64  a 

0.635 

3 

6 

10 

420 

64* 

0.649 

4 

8 

11 

462 

64  * 

0.672 

5 

4 

10 

420 

128* 

0.663 

6 

12 

10 

420 

64* 

0.692 

10 

10 

420 

128* 

0.764 

7 

13 

12 

504 

128* 

0.659 

13 

12 

504 

128* 

0.705 

8 

14 

9 

378 

128* 

0.652 

11 

8 

336 

128* 

0.546 

9 

'5 

10 

420 

128  * 

0.614 

10 

16 

9 

378 

256* 

0.699 

5 

9 

378 

128* 

0.587 

11 

14 

12 

504 

128* 

0.665 

11 

12 

504 

128* 

0.566 

We  then  computed  the  mean  P'(c)  for  each  Individual  subject  (on  the  basis  of  his  entire 
Series  B  performance)  The  results  are  listed  In  Table  IV  and  plotted  inthe  bar  graph  of  Figure 
6.  From  these  data  we  note  that  subjects'  performance  ranges  from  0.516  to  0  744.  In 
analyzing  this  variation  we  first  designated,  by  x  at  the  top  of  Figure  6,  those  listeners  who 
attended  nine  or  more  of  the  16  data- loop  tests  In  Series  B.  The  unmarked  subjects  attended 
only  one  to  four  of  the  tests.  We  suggest  that  the  results  of  the  eight  listeners  In  the  latter 
category  should  be  neglected  because  of  Insufficient  sampling. 
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The  listeners  were  instructed  to  listen  lor  certain  auditory  stimuli  and  judge  them  on  the 
basis  of  specified  criteria;  that  is,  they  were  to  use  a  weighted-criteria  decision  technique. 
Subjects  14M  and  17M  did  not  use  such  a  system,  and  we  have  evidence  that  subject  11F  (who 
attended  the  minimum  of  9  tests)  also  failed  to  apply  this  technique.  Oi  this  basis  we  recom¬ 
puted  the  mean  P'(c)  for  group  averages  of  Table  HI,  using  only  the  results  of  listeners  1  through 
8.  getting 

P'(c)  =  0.67“  t0.057  (13) 

Results  for  Series  A  and  Series  II  can  now  be  compared  on  the  basis  of  Equations  11  and  13 


TABLE  IV.  MEAN  P’(c)  FOR  EACH  SUBJECT  — SERIES  B 


Subject  P’(c) 


1M 

0.744 

±0.115 

2F 

0.709 

±  0.089 

3F 

0.701 

±0.117 

4M 

0.681 

±0.092 

5N1 

0.656 

±0.080 

6F 

0.648 

±0.071 

Subject  P'(c) 

8M  0.633  1  0.091 

9M  0.629  ±0.046 

10M  0.627  ±0.079 

11 F  0.622  ±0.087 

12M  0.622  ±0.055 

13F  0  616*0.189 


Subject _ P'(c) 

15M  0.590  ±0.030 

16N1  0.580  ±0.000 

17M  0.5C7  ±0.082 

18F  0.542  ±0.132 

19M  0.516  ±0.027 


FIGURE  6.  P'(C)  FOR  19  SUBJECTS  IN  SERIES  B  TESTS.  The  letters  M  and  F  denote  male  and  female 
subjects.  The  x  at  the  top  Identify  those  subjects  who  attended  nine  or  more  of  the  16  series  B  tests. 
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2.5.  DISCUSSION 

So  far  we  have  evaluated  the  auditory  performance  by  the  percentage  of  correct  Identifica¬ 
tions,  using  an  rnwolghted  mean  P(c)  of  alt  subjects  In  Series  A  and  a  weighted  mean  P'(c)  of 
eight  selected  subjects  In  Series  B.  We  can  analyze  the  data  for  additional  Information  by  ex¬ 
amining  the  rating  responses  of  Series  B  In  greater  detail.  Though  It  Is  difficult  to  establish 
likelihood  ratio  data  for  each  of  many  complex  stimuli,  we  can  break  down  the  net  effect  Into 
strict  (-2,  -2)  and  lax  (.1,  -l)  criteria.  This  has  the  effect  of  translating  a  cut  along  the  de¬ 
cision  axis,  as  Indicated  In  Figure  5,  and  should  generate  operating  characteristic  data  of  the 
type  discussed  earlier.  A  table  was  made  up  for  each  of  the  19  listeners  to  summarize  total 
responses  In  each  of  the  categories  *2,  *1.  -1,  -2  for  both  subsets,  EQ  and  EX.  From  these 
entries  Individual  probabilities  and  cumulative  probabilities  were  determined.  The  results  of 
this  statistical  calculation  are  shown  In  Figure  ",  In  which  Pj.q(A)  Is  plotted  on  the  ordinate 
and  Ptx<A)  on  the  abscissa  of  double-probability  paper.  The  positive  diagonal  of  these  normal- 
normal  graphs  Is  the  chance  line.  For  clarity,  lines  are  not  drawn  through  Individual  points, 
but  In  general  the  data  for  listeners  1  through  8  do  approximate  a  linear  relationship  with  unity 
slope.  (These  parameters  must  be  referenced  to  the  deviate  scale  to  the  right  and  top  of  the 
graph.)  However,  the  points  for  subjects  9  through  13  become  erratic,  and  those  for  subjects 
U  through  19  approach  the  chance  line.  We  are  remtnded  about  statements  made  earlier  con¬ 
cerning  listeners  9  through  19:  that  eight  of  them  had  an  Inadequate  auditory  sample  and  at 
least  two  of  the  other  three  did  not  use  a  weighted-criteria  decision  technique. 

A  negative  semi-diagonal  Is  drawn  on  the  graphs  of  Figure  "■  Oie  of  the  significances  of 
this  Is  that  the  Intersection  of  a  straight  line  (through  any  set  of  points)  with  this  diagonal  scales 
equal  probabilities  for  Identification  of  earthquakes  and  explosions.  Points  for  the  best  subjects 
were  connected;  the  resulting  ROC  curve  Is  plotted  on  linear  coordinates  In  Figure  8.  Observe 
from  this  curve  that  at  the  point  of  maximum  difference  between  hit  rate  and  false  alarm  rate 
the  percentage  of  correct  Identification  Is  about  67 -5X>  (equal  for  EQ  and  EX). 

2.6.  RESULTS 

In  the  Interpretation  of  these  data  the  question  whether  the  experiment  was  adequate  natur¬ 
ally  arises.  We  exposed  each  subject  to  an  average  of  1500  auditory  decisions  In  the  combined 
training  and  testing  program.  Is  this  sufficient?  It  Is  possible  that  the  rating  method  automates 
the  learning  process  so  that  the  listener  reaches  a  plateau  In  performance  more  rapidly  than  by 
the  two-choice  procedure  alone.  A  set  of  data  bearing  on  this  aspect  of  the  experiment  Is  pre¬ 
sented  In  Figure  9,  a  plot  of  the  results  throughout  the  Series  B  tests  for  11  listeners.  The 
points  are  relative  to  an  arbitrary  mean,  and  the  dashed  lines  are  least-square  lines  which  show 
how  each  subject’s  performance  varied  with  time  (time  and  experiment  are  synonymous  here). 


15 


Institute  of  Sdonco  and  Taehnotogy 


T  h •  University  af  Michigan 


MW  UAL  MVU1C 


M USUAL  MLUIt 


w»UH  KtUtl 


U»tUL  Mtutl 


FIGURE  7.  PIXIT  OK  Pkq(A)  vs  Pkx(A)  ton  19  SLIUKCTS.  PKq(A)  l»  plotted  against  PKX<A)  on  doul.le- 
probablllty  paper  The  slope  of  the  t>osi-rit  straight  Ur  •  to  Individual  sets  of  points  Is  determined  from  the 

deviate  scales. 


A  positive  slope,  such  as  that  for  observer  4M,  indicates  that  the  listener  was  stili  improving 
at  the  end  of  the  experiment.  A  negative  siope,  such  as  that  for  observer  3F,  means  the  op¬ 
posite.  We  note  that  the  results  for  most  of  the  subjects  are  flat  (in  the  least -squares  sense). 

Another  question  concerns  the  rather  wide  variation  between  performance  of  Individuals. 
Part  of  this  difference  is  obviously  due  to  the  Inadequate  sampling  of  some  subjects.  We  ob¬ 
tained  Malco  audiograms  of  13  listeners,  looking  for  possible  correlations  of  auditory  results 
with  frequency  variations  in  individual  hearing  levels.  These  data  are  shown  In  Figure  10. 
Most  of  the  audiograms  are  near  "normai."  One  of  the  better  subjects,  4M,  has  a  hearing  ioss 
of  about  40  db  per  octave  above  1000  cps.  No  correlation  is  evident. 
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KIGITtK  8.  HOC  CITtVE  I'Olt  THE  BEST  St'BJKCTS.  The  |>o»IUvc  illngonnl  Is  the 
chnnee  line.  The  slope  of  (he  curve  st  n  particular  operating  point  Is  related  to  the 
weighting  factor  u)  Introduced  In  Kquntlon  8 


An  interesting  observation  which  can  be  given  only  cursory  examination  in  this  report  is 
shown  in  the  data  of  Table  V.  Most  of  the  seismograms  used  in  the  study  were  vertical-compo¬ 
nent  ones,  but  this  tabic  contains  exceptional  cases  where  horizontal  components  were  also  time- 
compressed.  Column  P(c)  contains  listening  results  for  the  cases  when  the  channels  were  an¬ 
alyzed  separately.  Column  Pp(c)  contains  the  results  of  dual-component  listening  where  a  ver¬ 
tical  component  and  a  horizontal  component  were  simultaneously  played  through  separate  speak¬ 
ers.  In  the  few  cases  tested,  auditory  recognition  Is  Improved  by  the  addition  of  the  second 
dimension. 


Finally,  it  would  be  significant  to  compare  the  auditory  analysis  method  to  an  analytical 
technique,  both  applied  to  the  same  seismic  data.  Many  of  the  seismograms  used  In  this  report 
were  studied  by  others  [5,  6],  who  investigated  the  discrimination  of  seismic  signals  on  the 
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EXPERIMENT 

FIGURE  9.  CIIRONOUOGICAI,  RESULTS  FOR  11  SUBJECTS  IN  SERIES  R  TESTS. 
Points  arc  referenced  to  a  relative  mean,  ami  least-square  lines  (Hashed)  show  the 
variation  in  performance  with  experiment  or  with  time.  A  positive  slope  Indicates 
that  subjects  are  still  improving  at  the  end  of  the  experiment.  A  negative  slope 
means  the  opposite. 
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FIGURE  10  MAICO  A  t’lllOGR  A  MS  FOR  13  SUBJECTS  Solid  awl  dashed  lines  are  for  rlRht  awl  left  ears. 

respectively. 


TABLE  V.  DUAL-COMPONENT  VS.  SINGLE -COMPONENT  LISTENING 


Recording 

Event  Distance 
(km) 

M 

Component 

P(c) 

PD(0 

Castle  Cliffs,  Utah 

N-Coll. 

217 

3.0 

Z 

0.907  1 

O  QAO 

Castle  Cliffs,  Utah 

N-Coll. 

217 

3.0 

T 

0.945) 

u.you 

Uamona,  Calif. 

EX-Q 

150 

Z 

0.146} 

n  iq  a 

Ramona,  Calif. 

EX-Q 

150 

R 

0.228  f 

U.liTl 

Luera  Ranch,  N.  Mex. 

EQ 

315 

2.2 

Z 

0.669 } 

n  i*  a  o 

l.uera  Ranch,  N.  Mex. 

EQ 

315 

2.2 

T 

0.633) 

Lucra  Ranch,  N.  Mex. 

EQ 

315 

2.2 

Z 

0.665} 

H  740 

Luera  Ranch,  N.  Mex. 

EQ 

315 

2.2 

R 

0.523 

U .  1 13 

Quarter  Master  V.  P.,  Arlz. 

EX-N 

235 

2.6 

Z 

0.528} 

0  AO  7 

Quarter  Master  V.  P.,  Arlz. 

EX-N 

235 

2.6 

T 

0.484) 

basis  of  the  ratio  of  compressions!  wave  amplitude  to  shear-surface  wave  amplitude.  One  way 
of  comparing  results  would  lie  by  likelihood  ratios.  In  Table  VI  we  Illustrate  a  technique  for 
calculating  likelihood  ratio  data  from  other  criteria  by  using  probabilities.  The  first  column 
lists  the  different  stimuli  or  divisions  of  a  criterion.  The  second  and  third  columns  give  the 
number  of  times  each  stimulus  Is  observed  for  EQ  and  EX.  These  columns  are  summed  and  the 
Individual  probabilities  determined  as  Indicated  In  the  fourth  and  fifth  columns.  Likelihood  ratio 
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TABLE  VI.  LIKELIHOOD  RATIO  CALCULATION 


Criterion 

EQ 

EX 

peq*x* 

PEX(x) 

Kx) 

1 

nl 

mi 

Vlnt 

nij/Entj 

nj/Enj-i-nij/wntj 

2 

n2 

'"2 

V-"! 

m^Lntj 

n2/ln1-i-m2/lm1 

3 

n3 

m3 

n3/lni 

nij/Enij 

lystii+iyrm, 

4 

n4 

m4 

n4/lni 

ny-nij 

n4/-ni-i-n,4/Lmi 

In. 

i 

i 

f  (x)  is  then  the  ratio  ol  the  corresponding  probabilities.  Table  VII  presents  the  results  ol  carry¬ 
ing  oul  this  calculation  on  seismograms  used  in  this  report.  In  addition  to  the  amplitude  ratio 
mentioned  above,  we  have  also  established  the  values  ol  f (x)  lor  event  magnitude  and  range. 

These  are  readily  compared  with  the  auditory  data  in  the  table.  The  significance  ol  f(x)  is  in¬ 
dicated  by  the  extent  to  which  it  deviates  from  unity  (either  ♦  or  -). 

One  objective  of  the  auditory  experiments  was  to  test  the  listeners’  performance  as  a  func¬ 
tion  of  seismic  propagation  distance.  In  Table  VIII  we  present  these  results.  Note  that  there 
indeed  is  some  variation  in  P'(c)  with  distance,  although  it  must  be  admitted  that  65%  of  the 
seismic  data  available  in  this  study  were  In  the  eplcentral  distance  range  from  100  to  4)0  km. 
it  is  interesting  to  note  also  (In  Table  VIII)  the  percentages  (p*)  obtained  by  Witits,  ct  ai.  [5], 
from  the  same  seismic  data;  p*  Is  the  percent  of  earthquakes  which  have  a  shear/compressionai 
wave  amplitude  ratio  greater  than  the  corresponding  ratio  for  nuciear  explosions.  The  last  re¬ 
sult  In  the  P'(c)  column  of  Table  VIII  was  obtained  In  auditory  experiments  when  distance  was 
allowed  to  vary  randomly  over  a  broad  range  (100  to  3000  km)  in  a  given  iistening  test  and 
time  compression  was  held  constant.  The  reader  is  reminded  that  the  p*  coiumn  depicts  EQ 
identification,  whereas  the  P'(c)  coiumn  depicts  approximately  equal  identification  of  EOand  EX. 


TABLE  VII.  COMPARISON  OF  ANALYSIS  TECHNIQUES 


Auidtory 

f(x) 

S/P  Max. 

Kx) 

M 

Kx) 

R(km) 

Kx) 

+2 

2.04 

<0.63 

2.24 

<1 

.48 

<100 

1.46 

+1 

1.34 

0.63-1.25 

0.77 

1-2 

1.54 

100-200 

1.11 

-1 

0.75 

1.25-2.5 

0.64 

2-3 

0.97 

200-400 

0.58 

-2 

0.49 

2.5-3.75 

2.08 

3-4 

0.69 

400-1000 

1.29 

3.75-5.0 

0.77 

4-5 

1.59 

1000-3000 

0.61 

5.0-10.0 

3.16 

5-6 

1.07 

3000-7000 

2.97 
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TABLE  Vin.  DISCRIMINATION  PERCENTAGE  AS  A 
FUNCTION  OF  EPICENTRAL  DISTANCE 

P’(c)  Eptcenlral  Distance  p* 

(%)  (lent)  (<*) 


71 

35-100 

50-100 

85 

68 

100-200 

100-200 

59 

71 

200-400 

200-400 

52 

59 

400-700 

400-1000 

67 

63 

700-1700 

70 

3000-7000 

60 

100-3000 

‘From  Relerence  5. 

3 

CONCLUSIONS 

Seismic  signals  Irom  approximately  200  sources  (earthquakes  and  explosions)  were  time- 
compressed,  and  an  experimental  training  and  testing  program  was  conducted  with  the  resulting 
audio  sounds.  Experimental  parameters  were  controlled  to  cover  various  combinations  of  seis¬ 
mic  variables,  but  the  Inhomogeneity  characteristic  of  typical  field  seismograms  was  preserved 
and  the  signals  were  presented  to  the  listeners  in  random  order. 

There  arc  several  ways  of  analyzing  this  kind  of  data  In  selecting  an  approach,  we  were 
guided  in  part  by  the  multiplicity  and  complexity  of  auditory  stimuli  associated  with  seismo¬ 
grams.  One  method  measured  the  observer's  performance  on  the  lasts  of  percentage  of  cor¬ 
rect  identifications.  On  this  basis,  the  mean  score  for  all  listeners  Is  al>oul  651,,  and  the 
weighted  score  for  the  eight  subjects  who  most  closely  adhered  to  the  experimental  demands 
Is  about  68 

An  alternative  technique  classified  the  observer's  responses  In  a  simple  four-category 
rating  system  which,  iastcalty,  partitioned  the  decisions  according  to  "strict''  and  'lax  '  criteria. 
On  this  basis,  the  optimum  listener  Identified  about  67.51)  of  the  earthquakes  and  explosions. 

This  result  Is  based  on  a  priori  probabilities  of  EQ  and  EX  being  0.5.  The  added  significance 
of  the  alternative  technique  Is  that  an  ROC  curve  was  generated.  Techniques  are  available  In 
the  literature  (e.g.,  References  7  and  8)  which  make  It  possible,  by  using  ROC  data,  to  translate 
the  result  stated  above  into  a  domain  where  other  a  priori  probabilities  exist. 

In  this  experiment  each  subject  was  required  to  make  approximately  1500  auditory  decisions. 
One  might  speculate  on  the  significance  of  a  much  more  Intense  auditory  program.  Ckir  data  do 
Indicate  that  most  of  the  listeners  apparently  reached  a  plateau  In  mean  performance. 
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We  conclude,  then,  that  the  optimum  trained  observer  using  a  subjective  technique  of  aud¬ 
itory  analysis  Identified  approximately  two  thirds  of  the  seismic  signals  consisting  of  a  random, 
inhomogeneous  set  of  200  earthquake  and  explosion  seismograms.  Each  of  these  events  was 
recorded  at  one  geographic  location.  If  multiple  statistically  Independent  recordings  of  the 
same  event  are  available  (say  from  a  global  distribution  of  seismic  stations)  then  in  theory  the 
curve  of  Figure  8  can  be  raised  3  db  for  each  twofold  Increase  In  the  number  of  stations. 

In  a  practical  application  of  the  auditory  technique,  if  the  inhomogeneity  of  seismic  data 
were  limited  by  using  recordings  from  only  one  station,  then  It  is  conceivable  that  the  result 
of  this  experiment  could  be  Improved.  This  is  particularly  expected  If  discrete  "unknown”  sig¬ 
nals  are  tested  by  comparison  with  known  reference  signals,  both  EQ  and  EX,  peculiar  to  the 
azimuth  and  range  estimated  for  the  "unknown"  ones. 
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Appendix 

DESCRIPTIONS  OF  SPEEDED-UP  DATA  LOOPS  USED  IN  THE  TESTS 


Channel 

No. 

Event 

Date 

Recording 

Location 

Source 

Location 

Distance 

(km) 

Seis¬ 

mometer 

Magnitude 

(M) 

TAPE  LOOP  NO.  4 

1 

EQ 

7/28/60 

Kingman,  Arlz. 

276 

VB 

3.1 

2 

N-EX 

Kingman,  Arlz. 

NTS,  Nev. 

265 

VHS 

3.4 

3 

EQ 

8/24/59 

Shepherd,  Mont. 

Hebgen  Lake, 
Mont. 

250 

VW 

1.8 

4 

N-EX 

Kingman,  Arlz. 

NTS,  Nev. 

265 

VHS 

2.2 

5 

C-EX 

10/13/60 

Fletcher,  Nev. 

NTS,  Nev. 

2b5 

VW 

3.8 

6 

N-EX 

Quarter  Master 
View  Pt.,  Arlz. 

NTS,  Nev. 

235 

VW 

2.6 

7 

EQ 

8/24.  59 

Shepherd,  Mont. 

Hebgen  Lake, 
Mont. 

235 

VW 

2.0 

8 

N-EX 

Climax  Claims, 
Arlz. 

NTS,  Nev. 

226 

VW 

4.8 

9 

Q-EX 

7/13/62 

Flborn  Quarry, 
Mich. 

Humboldt,  Mich. 

214 

VHS 

2.7 

10 

N-EX 

Climax  Claims, 
Arlz. 

NTS,  Nev. 

220 

VW 

2.2 

11 

C-EX 

10/13/60 

St.  George,  Utah 

NTS,  Nev. 

220 

TB 

3  4 

12 

N-EX 

Climax  Claims, 
Arlz. 

NTS,  Nev. 

215 

VW 

2.8 

13 

C-EX 

10/13/60 

St.  George,  Utah 

NTS,  Nev. 

220 

LB 

3.7 

14 

N-EX 

Climax  Claims, 
Arlz. 

NTS,  Nev. 

220 

TW 

2.8 

EQ  —  Earthquake 
N-EX  — Nuclear  explosion 
C-EX  — Chemical  explosion 
Q-EX — Quarry  explosion 
U -EX  — Underwater  explosion 
B —  1  -Sec  Benloff 
HS  —  1/2-Sec  Hall-Sears  HS-10 
W  —  1-Sec  Wlllmore  Mark  I 
ET  —  1-Sec  Electrotech  EV-17 
V  —  Vertical 
T  —  Transverse 
L  —  Longitudinal 
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Channel 

No. 

Event 

Date 

Recording 

Location 

Source 

Location 

Distance 

(km) 

Seis¬ 

mometer 

Magnitude 

(M) 

TAPE  LOOP  NO.  5 

1 

EQ 

1/26/59 

Hollister,  Calif. 

Mt.  Hamilton, 

95 

VW 

2.7 

Calif. 

2 

N-EX 

12/10/61 

Hope,  N.  Mex. 

Carlsbad. 

100 

VW 

4.0 

N.  Mex. 

3 

EQ 

6/30/62 

Winslow,  Ariz. 

160 

VHS 

0.9 

4 

N-EX 

Castle  Cliffs, 

NTS,  Nev. 

187 

VW 

4.2 

Utah 

5 

EQ 

8/24/59 

Shepherd,  Mont. 

Hebgen  Lake, 

240 

VW 

2.9 

Mont. 

6 

N-EX 

Quarter  Master 

NTS.  Nev. 

235 

VW 

2.6 

View  Pt.,  Arlz. 

7 

EQ 

12/13/61 

Magdalena, 

315 

VW 

2.2 

N.  Nex. 

6 

N-EX 

Pica,  Ariz. 

NTS.  Nev. 

350 

VW 

3.1 

9 

EQ 

1/26/59 

Hollister,  Calif. 

Owens  Valley, 

432 

VW 

4.0 

Calif. 

10 

N-EX 

Hollister,  Calif. 

NTS.  Nev. 

450 

VW 

3.1 

11 

EQ 

7/24/60 

St.  George,  Utah 

Cloudcroft, 

638 

VW 

4.5 

N.  Mex. 

12 

N-EX 

Rifle,  Colo. 

NTS,  Nev. 

786 

VW 

4.1 

13 

EQ 

5/19/62 

U  of  M  Well. 

Acapulco. 

3000 

VHS 

Mich. 

Mexico 

M 

N-EX 

Fiborn  Quarry, 

NTS.  Nev. 

3000 

VHS 

3.35 

Mich. 

TAPE  LOOP  NO.  6 

1 

EQ 

7/30/62 

Neapnlts,  Crete 

402 

VHS 

5.1 

2 

U-EX 

7/1/61 

Coopers  Mill, 

Offshore  Maine 

150 

VW 

Maine 

3 

EQ 

7/28/62 

Neapolts,  Crete 

128 

VHS 

2.7 

4 

U-EX 

8/26/61 

DeBclI’s  Ranch, 

Off  San  Clemente 

123 

VW 

4.1 

Calif. 

Isle,  Calif. 

5 

EQ 

4/23/62 

Concepcion, 

Chile 

310 

VHS 

3.65 

6 

C-EX 

7/1/58 

Star  Lake, 

St.  Lawrence 

100 

VW 

3.35 

New  York 

Seaway 

7 

EQ 

1/26/59 

Hollister,  Calif. 

Mt.  Hamilton, 

118 

VW 

3.0 

Calif. 
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Channel  Recording  Source  Sels- 


No. 

Event 

Date 

Location 

Location 

Distance 

(km) 

mometer 

Magnitude 

(M) 

8 

EX 

8/3/61 

Scope  Lake, 
Alberta 

Suffleld,  Alberta 

100 

VHS 

2.2 

9 

EQ 

5/14/62 

Rincon,  Puerto 
Rico 

90 

VHS 

1.35 

10 

Q-EX 

8/6/59 

Hendricks 
Quarry,  Mich. 

Rogers  City 
Quarry,  Mich. 

156 

vw 

0.2 

11 

EQ 

4/23/62 

Concepcion, 

Chile 

115 

VHS 

2.65 

12 

Q-EX 

8/4/59 

Hendricks 
Quarry,  Mich. 

Rogers  City 
Quarry,  Mich. 

156 

VW 

0.8 

13 

EQ 

8/25/59 

McLeod,  Mont. 

Hebgen  Lake, 
Mont. 

110 

vw 

2.3 

14 

Q-EX 

7/12/61 

Hendricks 
Quarry,  Mich. 

Rogers  City 
Quarry,  Mich. 

156 

VB 

5.6 

TAPE  LOOP  NO.  7 

1 

EQ 

8/25/59 

McLeod,  Mont. 

Hebgen  Lake, 
Mont. 

110 

VW 

2.4 

2 

Q-EX 

6/29/62 

Ftborn  Quarry, 
Mich. 

Cedarville 
Quarry,  Mich. 

70.2 

VHS 

1.9 

3 

EQ 

8/25/59 

McLeod,  Mont. 

Hebgen  Lake, 

Mont . 

110 

VW 

2.0 

4 

U-EX 

8/20/62 

Barnett  Chapel, 

Ky. 

Off  Cape 
Girardeau,  Mo. 

100 

VHS 

0.8 

5 

EQ 

4/22/62 

Concepcion, 

Chile 

103 

VHS 

3  4 

6 

U-LX 

8/20/62 

Barnett  Chapel, 
Ky. 

Off  Cape 
Girardeau,  Mo. 

100 

VHS 

0.8 

7 

EQ 

7/29/62 

Neapolis,  Crete 

167 

VHS 

3.2 

8 

U-EX 

7/1/61 

Coopers  Mill, 
Maine 

Offshore  Maine 

146 

VW 

2.4 

9 

EQ 

4/23/62 

Concepcion, 

Chile 

310 

VHS 

3.65 

10 

U-EX 

7/1/61 

Coopers  Mill, 
Maine 

Offshore  Maine 

165 

VW 

2.4 

11 

EQ 

8/13/61 

Mt.  Laguna, 
Calif. 

123 

VW 

2.4 

12 

U-EX 

8/30/61 

DeBell's  Ranch, 
Call!. 

Off  San  Clemente 
Isle,  Calif. 

123 

VW 

2.8 

13 

EQ 

8/9/61 

Chief  Peak, 

92 

VW 

1.4 

Caitf. 
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Channel 

No. 

Event 

Dale 

Recording 

Location 

Source 

Location 

Distance 

(km) 

Seis¬ 

mometer 

Magnitude 

(M) 

14 

U-EX 

9/2/61 

DeBell’s  Ranch, 
Calif. 

Off  San  Clemente 
Isle,  Calif. 

123 

vw 

1.6 

TAPE  LOOP  NO.  8 

1 

EQ 

8/13/61 

Ml.  Laguna, 

Calif. 

90 

vw 

1.4 

2 

N-EX 

Radio  Crystal, 
Nev. 

NTS,  Nev. 

150 

vw 

2.8 

3 

EQ 

8/18/61 

Mt.  Laguna, 

Calif. 

94 

vw 

0.8 

4 

U-EX 

8/20/62 

Barnett  Chapel, 
Ky. 

Off  Cape 
Girardeau,  Mo. 

97.4 

vw 

1.2 

5 

EQ 

4/20/61 

Hollister,  Calif. 

85 

vw 

6 

EQ 

1/28/59 

Hoilister,  Calif. 

Mt.  Hamilton, 
Calif. 

95 

vw 

2.7 

7 

EQ 

4/16/61 

Hollister,  Calif. 

120 

vw 

8 

EQ 

8/23/59 

Red  Lodge, 

Mont . 

Hebgen  Lake, 
Mont. 

150 

vw 

3.6 

9 

EQ 

4/23/61 

Panoche,  Calif. 

150 

vw 

1.7 

10 

EQ 

8/23/59 

Red  Lodge, 

Mont. 

Hebgen  Lake, 
Mont. 

155 

vw 

2.2 

11 

EQ 

7/28/60 

St.  George,  Utah 

87 

vw 

12 

N-EX 

Climax  Claims, 
Arlz. 

NTS,  Nev. 

212 

VHS 

3.8 

13 

EQ 

3/1/62 

Castle  Cliffs, 

Utah 

78 

vw 

0.5 

14 

EQ 

7/31/62 

Neapolts,  Crete 

142 

VHS 

2.8 

TAPE  LOOP  NO.  9 

1 

EQ 

4/24/62 

Concepcion, 

Chile 

48 

VHS 

2.55 

2 

Q-EX 

8/17/62 

Townsend,  Mont. 

Trident  Quarry, 
Mont. 

80 

VHS 

2.7 

3 

EQ 

4/16/61 

Hollister,  Calif. 

45 

VW 

0.5 

4 

Q-EX 

6/29/62 

Ftborn  Quarry, 
Mich. 

Cedarvlile, 

Mich. 

70.2 

VHS 

2.1 

5 

EQ 

7/28/62 

Neapotls,  Crete 

76 

VHS 

2.7 

6 

Q-EX 

8/11/61 

Nett  Lake,  Minn. 

Morton  Quarry, 

70 

VW 

0.3 

Minn. 
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Channel 

No. 

Event 

Dale 

Recording 

Location 

Source 

Location 

Distance 

(km) 

Seis¬ 

mometer 

Magnitude 

(M) 

7 

EQ 

4/24/62 

Concepcion, 

Chile 

50 

VHS 

2.8 

6 

N-EX 

Corn  Creek, 

N'ev. 

NTS,  Nev. 

75 

VHS 

0.4 

9 

EQ 

7/7/62 

Climax  Claims, 
Arlz. 

60 

VHS 

0.4 

10 

Q-EX 

8/13/62 

Townsend,  Mont. 

Wolf  Creek, 

Mont. 

60 

VHS 

0.5 

11 

EQ 

6/9/61 

Mauna  Loa, 

Hawaii 

32 

vw 

2.0 

12 

Q-EX 

7/6/62 

Flborn  Quarry, 
Mich. 

Ccdarvllle, 

Mich. 

70.2 

VHS 

1.2 

13 

EQ 

5/12/62 

Rincon,  Puerto 
Rico 

45 

VHS 

1.15 

H 

Q-EX 

8/10/61 

Nett  Lake,  Minn. 

Pierce  Mine, 
Minn, 

81.6 

VW 

0.3 

TAPE  LOOP  NO. 

10 

1 

EQ 

4/16/61 

Hollister,  Calif. 

200 

VW 

2.1 

2 

Q-EX 

7/13/62 

Flborn  Quarry, 
Mich. 

Humboldt  Quarry 
Mich. 

213.9 

VHS 

2.7 

3 

EQ 

8/22/59 

Edgar,  Mont. 

Hebgen  Lake, 
Mont. 

200 

VW 

1.9 

4 

U-EX 

8/3/62 

Clayton,  N.  C. 

Offshore  North 
Carolina 

215.4 

vw 

2.3 

5 

EQ 

8/24/59 

Shepherd,  Mont. 

Hebgen  Lake, 
Mont. 

210 

vw 

2.3 

6 

U-EX 

8/18/62 

Waverly,  Tcnn. 

Off  Cape 
Girardeau,  Mo. 

200 

VHS 

1.8 

7 

EQ 

5/30/62 

Coconut  Pi., 
Amer.  Samoa 

200 

VW 

4.5 

8 

N-EX 

Climax  Claims, 
Arlz. 

NTS,  Nev. 

217 

vw 

3.4 

9 

EQ 

5/31/62 

Coconut  Pt., 
Amer.  Samoa 

200 

VHS 

3.6 

10 

N-EX 

Castle  Cliffs, 
Utah 

NTS,  Nev. 

190 

VW 

4.2 

11 

EQ 

8/24/59 

Shepherd,  Mont. 

Hebgen  Lake, 
Mont. 

225 

vw 

2.3 

12 

N-EX 

Climax  Claims, 
Arlz. 

NTS,  Nev. 

226 

vw 

4.8 
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Channel 

No. 

Event 

Date 

Recording 

Location 

Source 

Location  Distance 

(km) 

Seis¬ 

mometer 

Magnitude 

(M) 

13 

EQ 

8/24/59 

Shepherd,  Mont. 

Hebgen  Lake, 

Mont. 

210 

VW 

2.5 

14 

U-EX 

8/18/62 

Waverty,  Tenn. 

Off  Cape 

Girardeau,  Mo. 

200 

VHS 

1.8 

TAPE  LOOP  NO.  11 

1 

EQ 

5/31/62 

Coconut  Pt., 

Amer.  Samoa 

315 

VHS 

4.3 

2 

U  -EX 

7/5/61 

Hope,  Maine 

Offshore  Maine 

284.4 

VW 

2.4 

3 

EQ 

5/29/62 

Coconut  Pt., 

Amer.  Samoa 

390 

VHS 

4.2 

4 

U-EX 

7/2/61 

Coopers  Mill, 
Maine 

Offshore  Maine 

300 

VW 

2.4 

5 

EQ 

4/23/62 

Concepcion, 

Chile 

310 

VHS 

3.65 

C 

N-EX 

Pica,  Ariz. 

NTS.  Nev. 

322 

VW 

3.1 

7 

EQ 

4/24/62 

Concepcion, 

Chile 

280 

VHS 

3.1 

8 

N-EX 

Willow  Springs, 
Artz. 

NTS,  Nev. 

271 

VW 

2.6 

9 

EQ 

1/28/59 

Hollister,  Calif. 

432 

VW 

4.0 

10 

N-EX 

Winslow,  Artz. 

NTS,  Nev. 

554 

VHS 

3.1 

11 

EQ 

7/24/60 

St.  George,  Utah 

700 

VW 

3.9 

12 

U-EX 

7/1/61 

Coopers  Mill, 
Maine 

Offshore  Maine 

277.7 

VW 

2.2 

13 

EQ 

7/24/60 

Kingman,  Arlz. 

650 

VW 

3.2 

14 

N-EX 

Rifle,  Colo. 

NTS,  Nev. 

780 

VW 

4.8 

TAPE  LOOP  NO.  12 

1 

EQ 

8/18/61 

Ml.  Laguna, 

Calif. 

93 

VW 

1.3 

2 

Q-EX 

7/6/62 

Fiborn  Quarry, 
Mich. 

Rogers  City 
Quarry,  Mich. 

140 

VW 

4.0 

3 

EQ 

4/23/61 

Hoittster,  Calif. 

150 

VW 

1.2 

4 

Q-EX 

9/20/62 

Fiborn  Quarry, 
Mich. 

Humboldt  Quarry, 
Mich. 

200 

VB 

1.3 

5 

EQ 

7/28/60 

St.  George,  Utah 

87 

VW 

1.3 

6 

Q-EX 

7/6/62 

Fiborn  Quarry, 
Mich. 

Rogers  City 
Quarry,  Mich. 

140 

VHS 

4.0 
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Channel 

No. 

Event 

Date 

Recording 

Location 

Source 

Location 

Distance 

(km) 

Seis¬ 

mometer 

Magnitude 

(M) 

7 

EQ 

5/28/62 

Coconut  Pt.t 
Arner.  Samoa 

140 

VHS 

3.8 

8 

N-EX 

12/10/61 

Hope,  New 

Mexico 

Carlsband,  New 
Mexico 

100 

VW 

4.0 

9 

EQ 

7/31/62 

Neapolls,  Crete 

116 

VHS 

2.7 

10 

Q-EX 

8/19/60 

Hendricks 
Quarry,  Mich. 

Rogers  City 
Quarry,  Mich. 

156 

VW 

1.7 

11 

EQ 

4/23/62 

Concepclbn, 

Chile 

127 

VHS 

2.25 

12 

Q-EX 

7/11  '62 

Flborn  Quarry, 
Mich. 

Rogers  City 
Quarry,  Mich. 

141 

VHS 

2.1 

13 

EQ 

8/16/61 

Mt.  Laguna, 

Calif. 

144 

VW 

1.4 

14 

N-EX 

Radio  Crystal, 
Nev. 

NTS,  Nev. 

189 

VW 

2.6 

TAPE  LOOP  NO. 

13 

1 

N-EX 

Radio  Crystal, 
Nev. 

NTS,  Nev. 

189 

VW 

2.6 

2 

N-EX 

Radio  Crystal, 
Nev. 

NTS,  Nev. 

189 

TW 

2.6 

3 

N-EX 

Castie  Cliffs, 
Utah 

NTS.  Nev. 

217 

VW 

3.0 

4 

N-EX 

Castle  Cliffs, 
Utah 

NTS,  Nev. 

217 

TW 

3.0 

5 

Q-EX 

2/23/62 

Ramona,  Calif. 

Eagle  Mt.  Mine, 
Calif. 

150 

VW 

6 

Q-EX 

2/23/62 

Ramona,  Calif. 

Eagle  Mt.  Mine, 
Calif. 

150 

LW 

7 

EQ 

8/16/61 

Mt.  Laguna, 
Caiif. 

144 

VW 

1.4 

8 

EQ 

8/16/61 

Mt.  Laguna, 
Calif. 

144 

N-SW 

1.4 

9 

EQ 

12/13/61 

Luera  Ranch, 

N.  Mex. 

315 

VB 

2.2 

10 

EQ 

12/13/61 

Luera  Ranch, 

N.  Mex. 

315 

TB 

2.2 

11 

EQ 

12/13/61 

Luera  Ranch, 

N.  Mex. 

315 

VW 

2.2 

12 

EQ 

12/13/61 

Luera  Ranch, 

N.  Mex. 

315 

LB 

2  2 
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Channel 

No. 

Event 

Date 

Recording 

Location 

Source 

Location 

Distance 

(km) 

Seis¬ 

mometer 

Magnitude 

(M) 

13 

N-EX 

Quarter  Master 
View  Pt.,  Ariz. 

NTS,  Nev. 

235 

VB 

2.6 

14 

N-EX 

Quarter  Master 
View  Pt.,  Ariz. 

NTS,  Nev. 

235 

TB 

2.6 

TAPE  LOOP  NO.  i 

4 

1 

EQ 

7/28/60 

Kingman,  Ariz. 

276 

VB 

2.2 

2 

C-EX 

9/27/60 

Kanab,  Utah 

NTS,  Nev. 

315 

VW 

2.2 

3 

C-EX 

10/13/60 

Kanab,  Utah 

NTS,  Nev. 

287 

VW 

3.5 

4 

EQ 

9/27/60 

Kanab,  Utah 

275 

VW 

5 

EQ 

4/16/61 

Hollister,  Caiif. 

315 

VW 

6 

N-EX 

Climax  Claims, 
Ariz. 

NTS,  Nev. 

220 

V11S 

1.8 

7 

EQ 

7/27 / 60 

Kingman,  Ariz. 

220 

VB 

1.4 

8 

N-EX 

Climax  Claims, 
Ariz. 

NTS.  Nev. 

225 

VHS 

3.8 

9 

EQ 

12/13/61 

Luera  Ranch, 

N.  Mex. 

315 

VW 

2.2 

10 

N-EX 

Climax  Claims, 
Ariz. 

NTS,  Nev. 

225 

VHS 

3.4 

11 

EQ 

7/28  '60 

St.  George,  Utah 

320 

VW 

12 

N-EX 

Climax  Claims, 
Ariz. 

NTS,  Nev. 

225 

VHS 

3.7 

13 

EQ 

1/28/59 

Hoiiistcr,  Caiif. 

432 

VW 

4.0 

14 

N-EX 

Winslow,  Ariz. 

NTS,  Nev. 

562 

VW 

4.8 

TAPE  LOOP  NO. 

15 

1 

EQ 

7/2/63 

Shepherd,  Mont. 

Denver,  Coio. 

720 

VET 

4.2 

2 

L'-EX 

7/18/63 

Copper  Harbor, 
Mich. 

Lake  Superior 

255 

VHS 

2.6 

3 

EQ 

7/7/63 

Shepherd,  Mont. 

Central  Utah 

750 

VHS 

4.9 

4 

Q-EX 

7/12/61 

Hendricks 
Quarry,  Mich. 

Babbit  Quarry, 
Minn. 

521 

VB 

3.7 

5 

EQ 

8/31/63 

Chicken  Springs, 
Oregon 

Monterey  Cty., 
Caiif. 

700 

mis 

4.3 

6 

U-EX 

7/25/63 

Eagle  Lake, 

Minn. 

Lake  Superior 

400 

VHS 

1.6 

7 

EQ 

7/7/63 

Shepherd,  Mont. 

Off  Coast  of  N. 

1400 

VHS 

4.7 

Caiif. 
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Channel 

No. 

Event 

Date 

Recording 

Location 

Source 

Location 

Distance 

(km) 

Seis¬ 

mometer 

Magnitude 

(M) 

8 

N-EX 

Magdalena,  N. 

Mex. 

NTS,  Nev. 

840 

VB 

2.8 

9 

EQ 

8/30/62 

Marysville, 

Calif. 

Cache  Creek, 
Utah 

850 

SP-Z 

Lo 

5.5 

10 

N-EX 

Bonneville,  Wyo. 

NTS,  Nev. 

1000 

VHS 

3.6 

11 

EQ 

9/16/63 

100  Mile  House, 

B.  C. 

Kern  County, 
Calif. 

1760 

SP-Z 

Hi 

5.0 

12 

N-EX 

Idaho  Springs, 
Colo. 

NTS,  Nev. 

1100 

VHS 

4.8 

13 

EQ 

9/5/62 

Marysville, 

Calif. 

Cache  Creek, 
Utah 

850 

SP-Z 

Hi 

4.6 

14  N-EX 

TAPE  LOOP  NO.  16 

Truth  or  Conse¬ 
quences,  N.  Mex. 

NTS,  Nev. 

940 

VW 

2.8 

1 

EQ 

7/8/63 

Shepherd,  Mont. 

Mexico- 

Guatomala 

3625 

VHS 

4.5 

2 

N-EX 

Beverly,  Ohio 

NTS,  Nev. 

3000 

VHS 

5.5 

3 

EQ 

7/8/63 

Shepherd,  Mont. 

Kodiak  Is., 
Alaska 

4500 

VHS 

4.6 

4 

N-EX 

Fiborn  Quarry, 
Mich. 

NTS,  Nev. 

3000 

VHS 

5.5 

5 

EQ 

7/8/63 

Shepherd,  Mont. 

Colombia 

5550 

VHS 

4.0 

6 

N-EX 

U  of  M  Well, 

Ann  Arbor, 

Mich. 

NTS,  Nev. 

3000 

VHS 

5.5 

7 

EQ 

7/9/63 

Shepherd,  Mont. 

Costa  Rica, 
Panama 

4750 

VHS 

5.1 

8 

N-EX 

Fiborn  Quarry, 
Mich. 

NTS,  Nev. 

3000 

VHS 

4.7 

9 

EQ 

8/29/63 

Berry  Creek, 

Nev. 

Off  Coast  of 
Peru 

6400 

VHS 

6.5 

10 

N-EX 

Washington, 

Georgia 

NTS,  Nev. 

3000 

VHS 

5.5 

11 

EQ 

7/1/63 

Shepherd,  Mont. 

Kurile  Islands 

7000 

VHS 

4.5 

12 

N-EX 

Shepherd,  Mont. 

Fallon,  Nev. 

1100 

VHS 

4.9 

13 

EQ 

9/1/63 

Chicken  Springs, 
Oregon 

Nicaragua 

4500 

VHS 

4.4 

14 

N-EX 

Idaho  Springs, 
Colo. 

NTS,  Nev. 

1100 

VHS 

3.4 
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